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1  Overview 
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The  driving  idea  of  this  research  is  that  there  are  some  general  relationships  among 
dynamical  characteristics  of  vibrating  systems.  We  arrive  at  this  idea  when  we  com¬ 
pare  the  systems  studied  in  statistical  mechanics  and  structural  mechanics.  One  of 
the  favorite  systems  of  statistical  mechanics  is  a  set  of  mass  particles  connected  by 
springs  (Fig.  1),  We  obtain  the  same  kind  of  systems  if  we  discretize  the  equations  of 
elastic  structures.  For  example,  discretization  of  the  dynamical  equations  of  beams 
leads  to  a  chain  of  mass  particles  (Fig.  2),  while  for  more  complex  structures  one 
obtains  configurations  in  space  like  the  one  shown  in  Fig.  3. 

Another  example  of  the  similarity  is  black  body  radiation:  dynamical  equilibrium 
of  an  electromagnetic  field  in  a  closed  region  (symbolically  shown  in  Fig.  4).  Electro¬ 
magnetic  field  is  described  by  linear  Maxwell  equations.  It  interacts  with  the  charges 
on  the  wall  in  a  nonlinear  way.  There  is  a  very  close  structural  analogy  to  this  prob¬ 
lem:  a  linear  elastic  beam  connected  to  two  nonlinear  springs.  The  linear  beam  is 
described  by  linear  wave  equation  (which  mimics  Maxwell  equations). 

It  seems  natural  to  assume  that  the  laws  of  statistical  mechanics  might  be  valid 
for  structural  vibrations 

There  is,  however,  an  important  difference  between  the  ensembles  of  oscillators  in 
statistical  mechanics  and  structural  oscillators.  In  structures  we  often  have  to  take 
into  account  friction,  while  in  statistical  mechanics  oscillators  are  nondissipative. 
We  will  return  to  this  point  in  Section  6,  and  now  let  us  focus  only  on  the  elastic 
nondissipative  free  structural  vibrations. 

The  first  attempt  to  apply  statistical  mechanics  to  a  chain  of  oscillators  was  un¬ 
dertaken  about  40  years  ago  by  E.  Fermi,  J.  Pasta,  and  S.  Ulam  [1].  In  terms  of 
beam  vibrations,  they  considered  a  discretized  version  of  the  equation  for  nonlinear 
longitudinal  vibrations  of  an  elastic  beam  with  the  stress-strain  relation  of  the  form 


a  =  EY{e  +  pe^)  (1) 

where  a,  Eyi  £,  and  jS  are  stress.  Young’s  modulus,  strain  and  nonlinearity  parameter, 
respectively.  The  displacements  u{x,t)  obey  the  equation 


d'hi  _  ^ 


(2) 


where  i  is  time-(-^)2,  and  p  is  density  of  the  beam. 

It  was  a  common  belief  that  the  laws  of  statistical  mechanics  are  valid  for  systems 
with  a  very  large  number  of  degrees  of  freedom.  Since  a  1-D  continuum  possesses 
an  infinite  number  of  degrees  of  freedom,  one  may  expect  that  the  laws  of  statistical 
mechanics  can  be  applied  to  solutions  of  equation  (2).  To  study  the  behavior  of  these 
solutions,  equation  (2)  was  replaced  by  a  finite  difference  equation.  The  simplest 
test  for  validity  of  statistical  relations  is  the  equipartition  of  energy.  If  xii  is  the 
displacement  of  the  zth  node,  then  the  equipartition  of  energy  means 


0 


Figure  4:  (a)  Black  Body  Radiation,  (b)  Related  Structural  Problem;  Linear  Beam 
linking  Two  Nonlinear  Springs 


<  mul  >=<  mill  >=  ...  =<  mill  >  (3) 

where  m  is  the  node  mass,  and  <  •  >  denotes  the  time  average  along  a  trajectory: 
for  any  function  p(u,-,  h,) 


1 

<  p{ui,Ui  >=  lim  -  /  p{ui{t),Ui{t))dt 

0-*oo  C7  ^0 


(4) 


For  systems  considered  in  statistical  mechanics  (ergodic  systems)  motion  is  chaotic 
and  the  average  values  do  not  depend  on  the  initial  data.  The  common  value  (3)  is 
called,  by  definition,  temperature  T.  Equipartition  of  energy  is  a  necessary  condition 
for  the  laws  of  statistical  mechanics  to  be  true. 

The  discretized  equation  (2)  was  studied  numerically  [1]  for  the  number  of  points 
in  the  a:-direction,  N,  equal  to  64.  Surprisingly,  it  turned  out  that  equipartition  does 
not  hold  while  the  system  exhibits  a  recurrent  motion.  The  situation  was  cleared  up 
to  some  extent  by  KAM  theory  [2].  According  to  KAM  theory,  for  low  excitation 
energies,  the  beam  vibrates  approximately  as  a  linear  body,  and  there  is  no  equipar¬ 
tition  of  energy  over  the  degrees  of  freedom.  This  is  exactly  what  has  been  observed 
in  numerical  experiments  by  Fermi,  Pcista  and  Ulam,  where  excitation  energy  was 
relatively  small.  KAM  theory  tells  nothing  about  large  energy  excitation.  In  this 
case  only  numerical  results  have  been  obtained  so  far. 

Unfortunately,  in  the  papers  published  on  this  subject,  one  comes  across  many 
contradictory  statements.  Putting  aside  the  detailed  discussion  of  the  results  reported 
previously  (this  will  be  done  in  a  publication),  herein  the  picture  emerging  from  our 
numerical  simulations  is  described.  One  detail  of  this  picture-the  existence  of  upper 
energy  threshold-  can  be  found  in  previous  publications  [3]-[7],  but  the  entire  picture 
is  presented  here  for  the  first  time. 
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2  Longitudinal  Vibrations  of  Elastic  Beams 

At  present  the  only  way  to  establish  the  applicability  of  statistical  mechanics  to  struc¬ 
tural  vibrations  is  to  conduct  numerical  or  laboratory  experiments.  Computer  exper¬ 
iments  have  been  preferred  for  this  study  because  they  provide  very  easy  parametric 
control.  Since  our  goal  is  to  achieve  a  qualitative  understanding  of  the  situation,  we 
have  to  choose  a  structural  model  which,  on  one  hand,  is  as  simple  as  possible,  and 
on  the  other  hand,  captures  the  main  features  of  nonlinear  vibrations.  Longitudinal 
vibrations  seem  ideal  for  these  purposes. 

Longitudinal  vibrations  are  characterized  by  one  kinematical  field:  the  longitudi¬ 
nal  displacement  u{x,t).  Lagrangian  of  longitudinal  vibrations  has  the  form 

L  =  A{^pu^-Ui€))  (5) 

where  =  ^,  e  =  ^,  and  A  is  the  area  of  the  beam  cross-section.  Coordinate  x  is 
considered  as  the  Lagrangian  coordinate  of  beam  particles.  Energy  density  U{e)  is 
assumed  to  be  a  convex  function  of  strain  e.  The  corresponding  dynamical  equation 
of  free  vibrations  is 


d'^u  _  d  dU{Ux) 

^  dx  dux 

For  the  Fermi-Pasta-Ulam  (FPU)  problem 

U  =  '-Ey{s^  +  1&1)  (7) 

Expre.ssion  for  energy  density  (7)  can  only  be  used  for  moderate  strain  s,  because 
L  should  have  a  singularity  at  the  strain  value  £■  =  —!.  This  value  corresponds  to 
the  collapse  of  the  material  segment  to  a  point.  The  simplest  model  that  satisfies  the 
singularity  condition  is  the  Neo-Hookian  material 


/.L  2  3  -f-  Ux 


(8) 


where  p  is  the  shear  modulus  for  small  strains.  The  model  (8)  describes  the  elastic 
properties  of  some  rubbers  quite  well. 

We  performed  numerical  simulations  for  both  cases  (7)  and  (8).  They  look  quali¬ 
tatively  similar,  therefore  all  figures  presented  below  are  for  the  FPU  model.  The  only 
exception  is  Section  5,  where  we  evaluate  the  range  of  strains  where  our  conclusions 
are  valid  for  rubbers. 


To  set  up  a  finite  dimensional  model  we  discretize  the  continuum  beam,  and 
consider  a  chain  of  interacting  ‘‘particles””  (nodes)  with  the  Lagrange  function 


;V 


1 


=  E(2''‘ 

i=i  ^ 


U{ei))a 


(9) 


I 


where  u,-  is  the  displacement  of  the  ith  particle,  a  is  the  distance  between  neighboring 
particles,  and  -Ui/a,  a  -  N  =  £. 

We  use  periodic  boundary  conditions  for  our  simulations.  This  means  that  UAr+i  = 
tti.  In  fact,  our  results  do  not  depend  on  the  choice  of  the  boundary  conditions,  this 
will  be  explained  in  Section  5.  Discretized  version  of  the  equation  of  motion  is  taken 
in  the  form 


^ 

dui  dt  dui 


(10) 


The  laws  of  statistical  mechanics  were  derived  at  the  so  called  thermodynamical 
limit.  This  means  that  number  of  particles  and  energy  per  unit  volume  are  kept 
constant  while  the  size  of  the  system,  the  total  energy  and  the  total  number  of  particles 
N  tend  to  infinity.  In  our  case  this  corresponds  to  fixed  particle  distance  a  and  fixed 
energy  per  unit  length,  while  the  beam  length  £  =  aN  and  the  total  energy  E 


E  =  j‘(\puUU(e))dx  (11) 

tend  to  infinity. 

One  might  consider  another  limit,  continuum  limit,  when  the  length  of  the  beam  £ 
and  total  energy  E  are  fixed,  while  the  total  number  of  particles  N  tends  to  infinity.  In 
this  case  the  distance  between  neighboring  particles  a  =  £IN  tends  to  zero.  In  general, 
these  two  limits  are  different.  However,  longitudinal  beam  vibrations  possesses  the 
remarkable  property  that  the  thermodynamical  limit  coincides  with  the  continuum 
limit.  This  can  be  seen  if  we  scale  time,  axial  coordinate,  and  displacement,  and 
introduce  dimensionless  quantities 


where 


\W  can  write 


Ey 


=  lim 
£— 0 


U{£) 

ic2 

2- 


U{£)  =  EyUo{e) 


(12) 


where  Uo{e)  is  the  dimensionless  energy  density.  The  action  functional  of  the  system 
takes  the  form 


r  Ldt  =  r  AEY£LodT  Lo  =  -  Uo{v^)  (13) 

Jti  Jr\  Z 

Since  after  scaling  the  action  functional  has  the  same  form  (up  to  a  factor),  the 
equation  of  motion  also  has  the  same  form,  and  thermodynamical  limit  £  =  aN  oo 
is  reduced  after  scaling  to  continuum  limit  on  the  unit  segment  0  <  ^  <  1. 
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Note  that  the  coincidence  of  the  thermodynamical  and  continuum  limits  takes 
place  only  for  longitudinal  vibrations.  These  limits  are  different  for  lateral  vibrations. 
This  will  be  further  explained  in  Section  6. 

For  definiteness  the  following  will  be  carried  in  in  the  continuum  limit. 

One  useful  note  should  be  made  here.  In  the  case  of  the  FPU  model  the  scaled 
Lagrangian  has  the  form 


An  additional  scaling  is  possible  which  reduces  the  FPU  Lagrangian  to  a  standard 
form  with  /5  =  1 .  This  scaling  is 


we  have 


w 


T  — >  A  : 


Hence,  after  scaling,  dynamical  equations  of  the  FPU  model  coincides  with  the 
original  equation  with  /?  =  1.  Therefore,  all  results  are  presented  below  for  the  case 


/?=  1. 


3  A  Qualitative  Picture  of  Nonlinear  Vibrations 


3.1  Deviation  from  Equipartition 

Our  primary  goal  is  to  establish  whether  the  laws  of  statistical  mechanics  can  be 
applied  to  structural  vibrations.  To  this  end,  we  started  from  probes  of  equipartition 
of  energy  over  degrees  of  freedom.  We  measured  “partial  temperatures” 


9  1l 

Ti  =  pa  <  iq  >=  lim  -  /  paii^dt 
O^oo  U  Jo 

In  reality,  we  measured  the  quantities 


1 

0 


paufdt 


(14) 


(15) 


and  used  that  time  of  observation  6  can  be  taken  not  very  large  because  the  partial 
temperatures  approach  their  limit  values  fairly  fast.  Typical  dependence  of  quantities 
(15)  on  0  is  shown  in  Fig.  5.  In  all  the  graphs  time  refers  to  dimensionless  time  r 
(12).  The  largest  dimensionless  period  of  vibrations  of  a  linear  spring  is  of  order  unity. 
Usually  time  0  =  1000  is  enough  to  reach  the  limit  values  of  partial  temperatures. 
Deviation  from  equipartition  can  be  measured  by  the  quantity 
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A  =  max[ 


Ti-f 


Ti-f, 


T  T 

where  T  is  the  temperature  averaged  over  the  nodes 


T 


T1+T2  +  ...  +  T^v 

N 


3.2  Energy  Thresholds 

We  studied  the  dependence  of  A  on  the  energy  of  initial  excitation.  revealed  the 
existence  of  a  number  of  characteristic  values  of  energy,  energy  thresholds,  exceeding 
of  which  is  characterized  by  changes  in  beam  dynamics.  These  thresholds  are  shown 
in  Fig.  6.  Now  lets  discuss  what  all  the  lines  in  this  figure  mean. 

Lets  consider  an  initial  excitation  of  a  single  mode.  In  this  case  the  initial  data  is 
taken  in  the  form 


^t(O)  =  q-  sin 


f'i(O)  =  Pi  sin 


2'Kik 


N 

2'Kik 


N 


+  cos 
+  pf  cos 


2Kik 


N 

2Kik 


N 


(16) 


where  k  is  the  mode  number.  Hence  for  each  single  mode  of  excitation  we  have  4 
parameters  which  can  be  varied  arbitrarily  or  3  parameters  if  the  energy  value*  is 
prescribed. 

Lets  first  excite  the  lowest  mode.  If  the  energy  of  excitation  is  small,  the  final 
temperature  distribution  depends  significantly  on  initial  data.  'However,  if  we  exceed 
some  energy  value,  the  first  energy  threshold,  dependence  on  initial  data  disappears 
"while  energies  are  approximated  ecjually  distributed  over  the  nodes.  The  dependence 
of  A  on  the  energy  of  excitation  can  be  seen  in  Fig.  7.  There  is  a  fairly  sharp  drop 
iri  A  for  £  ~  10  .  Deviations  from  equipartition  A  for  E"  >  10~^  are  not  zero  but 

small,  in  the  order  of  0.05. 

In  fact,  we  observe  a  narrow  range  of  energies  for  w^hich  this  drop  occurs.  Moreover, 
if  we  try  many  initial  conditions  for  the  first  mode  excitation,  we  get  a  distribution 

^Energy  shown  in  all  the  graphs  is  dimensionless  energy  E", 


it  relates  to  the  “usual'’  energy  E  by  the  formula 


1 
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of  energy  thresholds  with  a  peak.  For  the  sake  of  brevity,  we  will  call  it  as  “energy 
threshold” ,  and  attach  the  corresponding  characteristic  value  of  energy. 

Simulation  show  that  the  value  of  the  energy  threshold  decays  with  the  growth 
in  the  number  of  degrees  of  freedom.  It  seems  very  plausible  that  the  first  energy 
threshold  reaches  some  asymptotic  value  w’hen  N  oo. 

What  happens  if  we  initially  excite  a  higher  mode,  say,  the  third  one?  The 
situation  is  very  similar  to  the  previous  one  with  the  only  difference:  energy  threshold 
for  the  third  mode  is  approximately  10  times  higher  than  the  first  energy  threshold. 
It  seems  that  it  also  approaches  some  limit  value. 

For  the  even  modes  the  behavior  of  the  energy  thresholds  is  different:  energy 
threshold  does  not  grow  monotonically  with  energy  growth.  We  assume  that  this 
degeneracy  relates  to  the  resonance  existing  between  even  modes  due  to  the  evenness 
of  the  number  of  degrees  of  freedom.  This  degeneracy  should  disappear  when  w'e 
consider  systems  with  an  odd  number  of  degrees  of  freedom.  This  hypothesis  has 
been  checked  for  17  degrees  of  freedom,  and  the  simulations  support  it  very  well. 

It  turns  out  that  there  are  two  energy  thresholds,  low  energy  threshold, and 
upper  energy  threshold,  such  that  for  energy  values  below  Ec  motion  is  ordered, 
and  for  values  above  E‘^  motion  is  approximately  ergodic,  and  equipartition  holds.  In 
the  intermediate  region  the  character  of  motion  depends  on  which  modes  are  excited. 
Low  energy  threshold  is  approximately  equal  to  the  first  energy  threshold. 

Low  energy  threshold  seems  not  to  depend  on  N  for  large  iV,  while  upper  energy 
threshold  seems  to  grow  with  increasing  N  although  for  small  N  we  observe  decay  of 
the  upper  energy  threshold  (see  Fig.  6). 

What  happens  if  we  excite,  for  example,  the  first  three  modes  simultaneously?  Of 
cour.se  the  energy  threshold  is  E3.  If  all  modes  are  involved  in  the  initial  excitation, 
then  equipartition  can  only  be  guaranteed  if  the  energy  of  excitation  is  greater  than 
the  upper  energy  threshold  E^^. 

Upper  energy  threshold  was  predicted  by  Chirikov  and  Israilev  [5]  and  was  ob¬ 
served  for  the  first  time  by  Bocchiery  et.  al.  [3]  for  a  chain  with  Lennard-. Tones 
interaction.  However  there  are  conflicting  statements  about  the  behavior  of  upper 
energy  threshold  when  N  grows.  For  example,  in  [5,Kantz]energy  thersholds  decay 
for  large  N.  As  it  was  mentioned,  our  results  support  simulations  by  Bocchiery 
et.al[3]. 


3.3  Energy  Transfer  from  Low  Frequency  Modes  to  High 
Frequency  Modes 

Let  us  now  consider  energy  thresholds  from  the  perspective  of  energy  transfer  between 
modes.  Let  only  the  first  mode  be  exited.  The  energy  spectra  for  various  values  of 
energy  are  shown  on  Fig.  8.  Do  not  pay  attention  to  the  “teeth”  on  the  curves,  as 
they  relate  to  the  short  time  of  averaging.  When  we  do  longer  runs,  these  curves  get 
smoother.  We  chose  a  time  of  averaging  in  such  a  way  that  we  obtain  results  within 
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a  reasonable  amount  of  time. 

For  low  energy  excitation,  E*  =  IQ-^,  which  is  far  below  energy  is  passed 
over  the  first  four  neighboring  modes.  When  E"  =  lO"”*,  first  8  modes  are  involved 
in  motion.  In  the  vicinity  of  Ec  (Fig.  8(e)),  about  half  of  the  modes  are  excited. 
This  seems  peculiar,  at  the  first  glance,  because  equipartition  of  spatial  distribution 
of  temperature  is  observed  at  these  values  of  energy.  It  is  known  from  statistical 
mechanics,  that,  if  equipartition  holds  in  one  set  of. coordinates,  it  also  holds  in  anv 
other  set  of  coordinates.  In  our  case  these  sets  of  coordinates  are  displacement  at 
nodes  and  amplitudes  of  modes.  The  explanation  relates  to  the  tail  of  the  curve  in 
Fig.  7.  Deviation  from  equipartition  is  small  but  not  zero.  This  “allows”  high  modes, 
in  Fig.  8,  not  to  be  excited.  Further  increase  of  energy  initiates  stronger  energy  flow, 
and,  when  energy  exceeds  the  upper  energy  threshold,  we  observe  equipartition  over 
modes  (Fig.  8(h)).  Note  a  specific  feature  of  energy  transfer:  decrease  of  energy  leads 
to  “freezing”  high  frequency  modes. 

3.4  Energy  Transfer  from  High  Frequency  Modes  to  Low 
Frequency  Modes 

Transfer  of  energy  from  high  modes  to  low  modes  (sometimes  referred  to  as  “backward 
energy  flow” ,  in  contrast  to  “forward  energy  flow”  in  the  case  of  energy  transfer  from 
low  to  higher  modes)  has  some  import  peculiarities.  To  discuss  them  consider,  for 
definiteness,  initial  excitation  of  the  5th  mode.  Energy  spectra  for  various  values  of 
excitation  energy^  are  shown  in  Fig.  9.  It  is  instructive  to  compare  Fig.  8  and  Fig.  9. 
First  of  all,  we  see  that  for  E'  =  10"^  only  the  fifth  mode  is  excited  (Fig.  9(a))  while 
the  first  mode  transmits  its  energy  to  the  three  neighboring  modes.  For  E^  =  10~^, 
the  resonance  of  fifth  and  tenth  modes  becomes  visible  and  only  these  two  modes  are 
excited  (Fig.  9(b)),  while  the  first  mode  conveys  its  energy  to  the  8  neighboring  modes. 
For  E  10  ,  the  resonance  of  fifth  and  tenth  modes  become  more  pronounced 

whde  the  energy  of  the  first  mode  is  transferred  to  the  12  neighboring  modes.  For 
=  10  ^  the  new  resonance  of  15th  mode  becomes  visible,  but  still  only  3  modes  are 
excited,  while  the  first  mode  shares  its  energy  with  15  other  modes.  We  conclude  that 
transfer  of  energy  occurs  very  easily  from  low  modes  to  high  modes,  and  backward 
energy  flow  is  impeded.  Physically,  it  seems  quite  natural  because  high  frequencv 
oscillators  are  more  stable  and  one  needs  to  supply  much  more  energy  in  order  to 
provide  the  interactions  with  other  oscillators.  Impediment  of  backward  energy  flow 
disappears  when  the  energy  exceeds  the  upper  energy  threshold;  energy  can  go  back 
and  forth  without  any  obstacles. 


^The  same  a.s  for  the  first  mode  on  Fig.  8. 


Figure  8;  Energy  Spectrum  for  First  Mode  Excitation 


Normalized  Correlation  of  a  Node  Velocity  and  the  Center  Node  Ve 
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3.5  Correlation  Radius 

In  the  preceding  subsections  we  have  seen  that  the  motion  is  not  completely  chaotic  at 
low  energy  threshold.  Another  confirmation  of  this  observation  comes  from  the  study 
of  the  correlation  radius.  We  measured  velocity  correlations  <  >  {i,j  =  1,2,. ..,N) 

and  considered  their  dependence  on  the  distance  between  the  nodes  i  —  j.  Typical 
graphs  of  this  dependence  are  shown  in  Fig.  10.  Here  <  iikUie  >  /  <  uie  >  is  plotted 
as  a  function  of  the  node  number  k  for  various  values  of  energy  of  excitation  of  the 
first  mode,  and  uig  is  the  displacement  of  the  center  of  the  32  particle  chain.  When 
E"  =  0.003,  which  is  below  Ec,  strong  correlation  is  observed.  Asymmetry  relates  to 
the  dependence  of  correlations  on  initial  conditions;  the  curve  retains  the  asymmetry 
of  the  prescribed  asymmetric  initial  data.  Around  the  threshold,  w’e  see  the  nonzero 
correlation  decays  as  energy  is  getting  larger.  When  the  upper  energy  threshold  E^  is 
reached  correlation  practically  disappears.  This  confirms  that  the  motion  for  E  >  E'^ 
is  practically  ergodic. 

3.6  Geometry  of  Phase  Space 

Above  discussion  sheds  some  light  on  the  geometrical  structure  of  the  phase  portrait 
of  beam  vibrations.  Unfortunately,  in  multidimensional  space  we  do  not  have  such 
powerful  tools  as  Poincare  maps  for  low-dimensional  systems.  Therefore,  we  have 
to  be  content  with  a  qualitative  picture.  The  sketches  trying  to  give  an  idea  on 
transformations  of  phase  portrait  are  presented  in  Fig.  11-13. 

If  the  energy  of  excitation  in  less  then  Ec,  then  the  energy  surface  is  occupied  by 
ordered  motion  (Fig.  11).  For  larger  energies  parts  of  the  islands  of  ordered  motion 
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decay  and  become  chaotic  “straits”  (Fig.  12).  Further  increase  of  the  energy  leads 
to  the  decay  of  the  “lowest  mode  island”  (Fig.  13).  Lowest  mode  island  disappears 
when  the  energy  exceeds  the  low  energy  threshold  Ec-  A  chaotic  sea  occupies  a 
considerable  part  of  the  energy  surface.  If  we  increase  the  energy  more,  the  chaotic 
sea  absorbs  the  second  mode  island  of  ordered  motion  and  so  on.  Of  course,  some  of 
the  islands  associated  with  modes,  can  disappear  simultaneously,  or  in  a  “inproper 
order”,  depending  on  the  peculiarities  of  the  system.  The  highest  mode  island  in  our 
simulations  survives  until  the  energy  exceeds  the  upper  energy  threshold.  Then  the 
chaotic  sea  occupies  the  entire  energy  surface,  except  very  small  islands  of  ordered 
motion,  and  motion  becomes  practically  chaotic. 

4  Is  There  a  “Planck  Constant”  for  Vibrations  of 
Beams? 

Consider  again  the  spectral  distribution  for  energy  of  excitation  lying  within  the 
interval  [Ec,  E^].  If  motion  were  ergodic,  energy  is  equally  distributed  over  the  modes, 
and  the  energy  spectrum  is  the  straight  horizontal  line  (Fig.  14(a)).  The  ordinate  of 
this  line  is  T,  because  all  modes  carry  the  energy  T  due  to  equipartition. 

In  reality,  we  observe  spectrum  distribution  shown  in  Fig.  14(b).  This  indicates 
that  motion  is  not  ergodic.  We  see  that  high  modes  are  practically  not  excited,  they 
are  ’’frozen.”  If  we  decrease  energy  further,  more  of  the  modes  become  frozen.  This  is 
very  reminiscent  of  the  difficulties  which  appeared  in  theoretical  predictions  of  heat 
capacity  of  gases  and  black  body  radiation  at  low  and  high  temperatures  100  years 
ago. 

A  theoretical  formula  for  a  spectrum  distribution  based  on  ergodic  hypothesis 
and  equipartition  was  suggested  by  Rayleigh  and  Jeans.  In  the  ID  case  it  is  exactly 
a  straight  line  as  shown  in  Fig.  14(a).  In  fact,  Rayleigh  and  Jeans  considered  the 
3D  case,  where  spectrum  is  parabolic  because  the  number  of  modes  in  an  interval 
+  du;]  is  proportional  to  u)^ .  Experiment  confirmed  the  Rayleigh- Jeans  formula 
for  low  frequencies,  but  revealed  strong  deviations  for  high  frequencies.  Qualitatively, 
the  Rayleigh-Jeans  spectrum  and  experimental  spectrum  are  shown  in  Fig.  15. 

Thinking  about  this  contradiction,  Jeans  noticed  that  high  frequency  motion  be¬ 
haves  as  if  it  were  frozen  when  we  lower  the  temperature  (or  energy). 

In  1900,  Planck  found  a  simple  formula  for  spectrum  distribution  which  fits  very 
well  with  the  experimental  data.  If  we  rewrite  this  formula  for  ID  vibration,  it  takes 
the  form 


E(u) 


hu) 


hui 

e  ^ 


where  h  is  a  constant  determined  from  fitting  to  experimental  data. 


(17) 
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For  low  modes,  T  >  hco,  the  exponent  in  (17)  can  be  expended  in  a  Taylor  series, 
and  we  obtain  equipartition 

E{u;)  =  T 

For  high  modes,  hui  ';$>  T,  we  have  exponential  decay  of  spectrum 

Eiu)  =  huje-^ 

The  form  of  Planck’s  spectrum  is  shown  qualitatively  in  Fig.  16. 

Later  on,  Planck  invented  the  ergodic  statistics  of  oscillators  which  leads  to  the 
energy  spectrum  (17).  The  statistics  were  based  on  quantum  hypothesis,  w'hich  in 
turn  gave  rise  to  quantum  mechanics. 

Comparing  the  curves  of  Fig.  8  and  Fig.  16,  we  see  a  remarkable  qualitative 
coincidence.  Moreover,  the  underlying  physics  seem  familiar;  high  frequency  modes 
freeze  if  we  lower  energy.  This  suggests  a  number  of  ideas: 

1.  Try  to  fit  experimental  data  for  beam  vibrations  using  Planck’s  distribution  and 
check  w'hether  the  Planck  constant  is  a  universal  one  for  vibrations  of  beams. 

2.  Check  to  what  extent  Planck’s  statistics  is  applicable  to  beam  vibrations. 

3.  We  know  that  the  dynamical  origin  of  the  spectrum  decay  is  nonergodicity  of 
motion.  Maybe,  there  is  an  explanation  of  behavior  of  the  quantum  oscillator 
based  on  classical  mechanics  and  nonergodicity  without  the  quantum  hypothe¬ 
sis? 

4.  Similarity  of  beam  spectrum  and  quantum  oscillator  spectrum  suggests  that 
some  correspondence  might  exist  between  statistical  mechanics  of  nonergodic 
systems  and  quantum  statistics. 

At  present,  we  are  far  from  understanding  all  of  that.  Trying  to  get  some  insight, 
we  started  from  an  attempt  to  fit  experimental  data  for  a  beam  energy  spectrum 
using  Planck’s  formula  (17).  A  typical  fit  is  shown  in  Fig.  17.  The  Planck  constant 
for  beam  vibrations  can  be  presented  in  the  form 

I),  =  Mclh  (18) 

where  M ,  c,  and  h  are  the  mass  of  beam,  speed  of  sound,  and  a  dimensionless  constant, 
respectively.  The  approximate  value  of  h  is 

h  =  0.003  (19) 

We  are  not  in  a  position  to  claim  that  h  is  a  universal  constant,  universal  in  the 
sense  that  it  is  only  material  dependent.  To  establish  that,  we  have  to  conduct  a 
two-parametric  study  with  respect  to  energy  of  excitation  and  the  number  of  degrees 
of  freedom.  W'e  have  not  done  this  yet  because  the  idea  of  using  Planck’s  spectrum 
came  only  recently.  At  present,  we  may  say  with  some  confidence  that  the  value  (19) 
does  not  depend  on  the  number  of  degrees  of  freedom.  A  comprehensive  study  is  now 
underwav. 
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5  High  Energy  Vibrations  of  a  Cantilever  Beam 

Now  we  are  going  to  give  some  additional  support  to  the  idea  that  statistical  mechan¬ 
ics  is  sensible  for  structural  vibrations. 

Consider  free  vibrations  of  a  cantilever  beam.  Let  us  excite  the  lowest  mode. 
Which  temperature  distribution  shall  we  get?  Temperature  at  the  clamped  edge 
should  be  zero;  temperature  at  the  free  edge  should  be  maximal.  Our  intuition, 
developed  on  linear  vibrations,  suggests  the  curve  shown  in  Fig.  18. 

But  this  is  in  obvious  contradiction  to  equipartition!  Therefore,  either  our  intu¬ 
ition  fails,  or  statistical  mechanics  does  not  work.  What  do  the  experiments  show? 

If  energy  of  excitation  is  small,  we  observe  a  temperature  distribution  which  really 
looks  like  the  curve  in  Fig.  18.  For  E  =  TO"®  it  is  shown  in  Fig.  19.  Note  that  motion 
is  not  ergodic  for  this  value  of  energy.  If  we  increase  the  energy  level  then  temperature 
distribution  changes;  for  E  =  10“^  it  is  shown  in  Fig.  20.  If  we  keep  increasing  energy 
and  overcome  the  low  energy  threshold,  temperatures  become  constant  and  we  get 
practically  ideal  equipartition  (Fig.  21). 

It  can  be  seen  in  Fig.  22  how  temperature  distribution  approaches  to  equilibrium 
distribution;  three  curves  correspond  to  successive  moments  of  observation. 

Again,  for  a  cantilever  beam  we  found  the  same  energy  thresholds  as  above:  if  the 
first  mode  is  excited,  transition  to  equipartition  occurs  at  the  first  energy  threshold; 
if  the  second  mode  is  excited,  then  for  equipartition  energy  should  exceed  the  second 
energy  threshold,  and  so  on.  Final  temperature  distribution  does  not  depend  on 
initial  excitation  ii  E  >  E'^. 

Another  conclusion  which  can  be  drawn  out  of  these  simulations  is  the  indepen¬ 
dence  of  the  equilibrium  state  from  boundary  conditions.  It  is  seen  that  the  edge 
zone  does  not  penetrate  into  the  interior  portion  of  the  beam.  This  is  why  our  results 
for  periodic  boundary  conditions  are,  in  fact,  universal. 

The  simulations  have  been  conducted  for  both  FPU  and  Neo-Hookean  models. 
For  Neo-Hookean  material  w'e  found  that  the  low'  energy  threshold  corresponds  to  the 
strains  of  order  0.2  —  0.3.  This  is  a  working  regime  for  many  rubbers  and  polymer 
materials. 


6  Dissipative  Systems 

Thermodynamics  and  statistical  mechanics  of  dissipative  system  of  oscillators  is  an 
open  field  now,  and  only  a  few  statements  (like  the  existence  of  attractors)  can  be 
made.  An  attempt  to  extend  thermodynamics  to  limit  cycles  was  undertaken  in 
papers  [12]  [13]  [14]  .  The  major  re, suits  are: 

1.  Entropy  disappears  from  thermodynamical  relations  because,  after  the  transi¬ 
tional  period,  vibrations  do  not  depend  on  initial  data.(initial  data  enters  in 
thermodynamical  relations  by  means  of  entropy 
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2.  The  averaged  Lagrangian  takes  place  of  energy  in  thermodynamical  relations. 

The  general  theory  [10]  has  been  considered  for  the  examples  of  Dufl&ng’s  Oscillator 

[13]  and  cantilever  beam  [14]. 

7  Major  Accomplishments 

1.  It  is  established  that  there  exists  a  sequence  of  energy  thresholds  for  beam  vi¬ 
brations.  They  characterize  the  energy  levels  which  should  be  exceeded  in  order 
to  get  equipartition  of  kinetic  energy  over  the  beam  particles.  The  correspond¬ 
ing  energy  level  depends  on  the  wave  length  of  initial  excitation.  The  lowest 
energy  threshold  Ec  corresponds  to  the  first  mode  excitation.  Second  energy 
threshold  E2  corresponds  to  the  second  mode  excitation,  and  so  on.  For  a  beam 
model  with  N  particles  there  are  y  energy  thresholds.  The  maximum  energy 
threshold  (upper  energy  threshold  E’^)  corresponds  to  the  transition  to  devel¬ 
oped  chaos.  For  E  >  E^  the  beam  motion  can  be  described  by  the  relations 
of  statistical  mechanics.  For  E  <  Ec  the  beam  motion  is  mostly  ordered.  In 
the  intermediate  region  Ec  <  E  <  E^  motion  of  the  beam  is  “partially  chaotic” 
in  the  sense  that  energy  exchange  between  beam  particles  is  enough  to  reach 
equipartition  but  not  enough  to  develop  a  completely  ergodic  motion. 

2.  The  existence  of  upper  energy  threshold  has  been  established  previously  in  a 
number  of  publications,  however,  there  are  contradictory  statements  on  the 
dependence  of  upper  energy  threshold  on  the  nummber  of  degrees  of  freedom. 
Our  simulations  show  that  upper  energy  threshold  decays  at  low  number  of 
degrees  of  freedom  and  seems  to  grow  at  large  number  of  degrees  of  freedom. 

3.  The  energy  transfer  has  been  studied  in  the  energy  range  Ec  >  E  >  E'^.  It 
was  observed  that  energy  is  easily  transferred  from  the  low  modes  to  higher 
modes  but  the  backward  energy  flow'  (from  high  modes  to  low  ones)  is  impeded. 
Analtically  the  energy  transfer  might  be  characterized  by  the  energy  spectrum 
w’hicli  seems  to  be  very  w'ell  fitted  by  the  Planck  spectrum. 

4.  In  computer  simulations,  the  approach  to  the  state  with  homogeneuous  distri¬ 
bution  of  temperature  has  been  clearly  observed  for  vibrations  of  a  cantilever 
beam  at  high  energies. 

5.  Thermodynamics  of  limit  cycles  has  been  developed  and  applied  to  Buffing’s 
oscillator  and  cantilever  beam 
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low  mode 

Figure  11;  Sketch  of  energy  surface  for  excitation  less  than  E^. 
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Figure  12:  Sketch  of  energy  surface  for  excitation  less  than  E, 
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'igure  21:  Temperature  Distribution  Along  the  Beam  with  E  = 
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Figure  22:  Evolution  of  Equilibrium  for  E  =  0.1 


0.1 


